Structural phase transitions in geometrically frustrated antiferromagnets 



00 

o 
o 

(N 
C 

o 



CO 



I 

C 

O 

o 



> 

m 

On 

in 
rn 

rn 
O 
00 

O 



T. E. Saunders and J. T. Chalker 
Theoretical Physics, University of Oxford, 1 Keble Road, Oxford, 0X1 3NP, United Kingdom 

(Dated: June 30, 2008) 

We study geometrically frustrated antiferromagnets with magnetoelastic coupling. Frustration 
in these systems may be relieved by a structural transition to a low temperature phase with re- 
duced lattice symmetry. We examine the statistical mechanics of this transition and the effects 
on it of quenched disorder, using Monte Carlo simulations of the classical Heisenberg model on the 
pyrochlore lattice with coupling to uniform lattice distortions. The model has a transition between a 
cubic, paramagnetic high-temperature phase and a tetragonal, Neel ordered low-temperature phase. 
It does not support the spin-Peierls phase, which is predicted as an additional possibility within 
Landau theory, and the transition is first-order for reasons unconnected with the symmetry analysis 
of Landau theory. Quenched disorder stabilises the cubic phase, and we find a phase diagram as a 
function of temperature and disorder strength similar to that observed in Zni-^Cd^C^CU. 

PACS numbers: 75.10.Hk 75.10.Nr 75.50.Lk 



I. INTRODUCTION 

Much of the interesting behaviour observed in geo- 
metrically frustrated antiferromagnets^ can be related to 
macroscopic ground state degeneracy in classical models 
for these systems. Some of the simplest models, with 
only nearest neighbour exchange interactions, remain in 
the paramagnetic phase down to zero temperature i 2 ' 3 ' 4 
Characterising the strength of exchange interactions by 
the magnitude of the Curie- Weiss constant ©cw, the 
paramagnetic phase is strongly correlated at tempera- 
tures T <C |@cw|- More realistically, ground state de- 
generacy is lifted by additional interactions or coupling 
between spins and other degrees of freedom. As a re- 
sult, the magnet may order or freeze at a temperature 
T c . If the relevant energy scales are small compared to 
exchange, T c <C |©cw|- O ne such escape route from the 
strongly correlated paramagnet arises from magnetoelas- 
tic coupling, which may lead to a magnetically driven 
cooperative Jahn- Teller transitions^ An alternative 
is that weak quenched disorder leads to spin freezing 
at low temperature^ In this paper we study structural 
phase transitions induced by magnetoelastic coupling in 
geometrically frustrated antiferromagnets and their sup- 
pression by quenched disorder. 

Recent theoretical workS£ on magnetic Jahn- Teller 
transitions in geometrically frustrated magnets has com- 
bined an exact treatment of a tetrahedral cluster of four 
spins with Landau theory applied to the pyrochlore lat- 
tice. Results for the cluster show that arbitrarily weak 
magnetoelastic coupling generates an instability to a dis- 
torted ground state. Landau theory elevates this instabil- 
ity to an ordering transition, with several possible phase 
diagrams!^ Depending on the distortion mode, the tran- 
sition from the cubic, paramagnetic phase to a tetragonal 
phase may be a continuous one, or it may be first order. 
If the transition is continuous, then spin Peierls order is 
expected in the tetragonal phase immediately below the 
critical point, but may give way to Neel order at lower 
temperature. Alternatively, if the cubic to tetragonal 



transition is discontinuous, it may either be to a phase 
with spin Peierls order, or directly to one with Neel order. 

While Landau theory provides an exhaustive listing 
of possibilities, a microscopic approach is necessary to 
decide which occur in a given model. Moreover, it may 
happen that transitions which on symmetry grounds need 
not be first order, because there is no cubic invariant in 
the Landau expansion for the free energy, in fact are dis- 
continuous, because a fourth-order invariant has a nega- 
tive coefficient. These considerations motivate the work 
we describe here. 

We examine as a representative model the classical 
Heisenberg antiferromagnet with nearest neighbour ex- 
change interactions on the pyrochlore lattice, includ- 
ing magnetoelastic coupling to distortions from cubic to 
tetragonal symmetry. Ordering occurs if a lattice distor- 
tion is accompanied by a free energy gain in the spin sys- 
tem that outweighs its cost in elastic energy. The order- 
ing temperature T c depends on the strength of the mag- 
netoelastic coupling compared to the appropriate elas- 
tic modulus. If coupling is weak, T c <C |@cw|- In the 
weak-coupling regime, thermodynamic behaviour is in- 
dependent of the strength of exchange interactions and 
a universal function of the ratio T/T c . We investigate 
the balance between magnetic free energy and elastic en- 
ergy in this regime by using Monte Carlo simulations 
of the spin system to provide a numerically exact treat- 
ment of its statistical mechanics. We combine this with a 
mean field treatment of the lattice degrees of freedom, in 
which we allow only uniform lattice distortions. We be- 
lieve our approach brings the advantages of a microscopic 
treatment of the magnetic degrees of freedom that drive 
ordering, without the complications and additional pa- 
rameters that would arise in a full treatment of phonons. 

Our results establish a rather different picture of the 
magnetic Jahn- Teller transition from that suggested pre- 
viously. We find that the structural transition is strongly 
first order, and to a Neel ordered state, independently of 
any cubic invariant in the elastic energy. In this sense, 
the transition is driven by the energy gain arising from 
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Neel order, and not from spin Peierls order. Since the 
jump we observe in the Neel order parameter is large, we 
believe this conclusion is likely to be robust, even though 
a change in microscopic details of the model may in prin- 
ciple favour spin Peierls order. 

The fact that arbitrarily weak magnetoelastic coupling 
is sufficient to generate a ground state instability prompts 
one to ask why structural phase transitions are not ubiq- 
uitous in geometrically frustrated antiferromagnets. We 
suggest here that weak exchange disorder (which is gen- 
erated by random strains in the presence of magnetoelas- 
tic coupling) acts to stabilise the high symmetry lattice 
structure. As we have shown elsewhere, in the absence 
of a structural transition, there is spin glass ordering. 9 In 
this sense, the magnetic Jahn- Teller and spin glass tran- 
sitions are alternative low-temperature fates for a geo- 
metrically frustrated antiferromagnet. 

The remainder of the paper is organised as follows. In 
Sec. m we define our model and summarise our main 
conclusions. We present a Monte Carlo study of the 
thermally-driven transition in a clean system in Sec. IIII1 
and consider the disorder-driven transition at zero tem- 
perature in Sec. HVl Finally, in Sec. [V] we discuss our 
results in relation to experiments on Zni_ K Cd x Cr204. 



II. MODEL AND SUMMARY OF RESULTS 

The model we treat has degrees of freedom of two 
types: classical Heisenberg spins Si of unit magnitude at 
the sites i of a pyrochlorc lattice; and a uniform tetrag- 
onal distortion of this lattice, parameterised by a single 
variable A. We include four contributions to the energy 
of the system: the exchange energy Ho of the undistorted 
lattice; the change in exchange energy Ha due to magne- 
toelastic coupling; the elastic energy H e \; and the energy 
Hdis arising from quenched exchange disorder. 

With interaction strength J the exchange energy of the 
cubic pyrochlore lattice is 

H = J^S i -S j , (1) 

ij 

where the sum is over nearest neighbour pairs. This 
Hamiltonian by itself is known to have a macroscopically 
degenerate ground state and to remain in the paramag- 
netic phase down to T = Our concern here is to find 
how other interactions lift this degeneracy and induce low 
temperature phase transitions. 

We express the change in exchange energy due to mag- 
netoelastic coupling in terms of the variable A. For def- 
initeness, we take the lattice distortion to be a compres- 
sion of the cubic z axis and an expansion of the x and 
y axes. We characterise the amplitude of this distortion 
by its effect on exchange interactions, taking the inter- 
action strength in the distorted system to be J — A be- 
tween neighbouring spin pairs in the same x-y plane, and 
J + A/2 between pairs in adjacent x-y planes. Then the 



magnetoelastic coupling is 

Wa X A <A- S - ( 2 ) 

where A,j has the values — A and A/2 for pairs i,j of the 
two types. The ground states of Ho + Ha for A > are 
Neel ordered: all spins in each x-y plane of the lattice are 
ferromagnetically aligned, with opposite orientations on 
successive planes, as illustrated in Fig. [T] (The opposite 
choice of sign, A < 0, leads to classical ground states 
with high degeneracy, and we do not consider it further.) 




FIG. 1: Classical groundstate spin configurations on a sin- 
gle tetrahedron. Left: for an undistorted tetrahedron all ex- 
change interactions are of equal strength and ground states 
include non-collinear spin configurations as illustrated. Right: 
for a tetragonally distorted tetrahedron, two longer bonds 
(broad lines) and four shorter bonds result in modulated in- 
teraction strengths, and ground state spin configurations are 
collinear. 

We include elastic energy within a harmonic approxi- 
mation, writing for a system of N spins 

H cl = ^NA 2 . (3) 

The coupling constant k is related to an elastic coeffi- 
cient c and the derivative dJ/dx of exchange interaction 
strength with distance in the following way. First, A can 
be expressed in terms of a strain e and an equilibrium 
bond length a, as the product A = eadJ/dx. Taking 
the elastic energy density to be ce 2 /2, we obtain Eq. ([3|) 
with k ~ ca/(dJ/dx) 2 . Thus k is large if magnetoelastic 
coupling is weak. 

Exchange disorder is represented by 

^dis = ^ SijSi ■ Sj (4) 

ij 

with Sij non-zero only for nearest neighbours, and drawn 
independently for each such pair i,j from a distribution 
that is uniform over the range [—8, 8]. 

Our objective is to study the model defined by the 
Hamiltonian 

H = Ho + Ha + Hd + Hdis (5) 
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in the regime n J^ 1 and A, 8 <C J. In these cir- 
cumstances Hq makes the dominant contribution to the 
energy of the system and transitions from the paramag- 
netic phase take place at temperatures T <C J (we set 
&b = 1 throughout this paper). In Fig. [2] we show a 
schematic phase diagram for the model as a function of 
the dimensionless variables kT and kS, which we now dis- 
cuss. First, consider the system without quenched dis- 
order (5 — 0) and at T = 0. Ground states are Neel 
ordered for A > and in these states Ha = —2NA. As 
a result, H is minimised by a non-zero tetragonal lattice 
distortion, with A = 2/k, no matter how weak the mag- 
netoelastic coupling and hence how large the coefficient k. 
The net energy gain per spin from this lattice distortion is 
2/k, and so one expects a transition to the paramagnetic 
phase at a temperature T c ~ kt x . As we demonstrate 
in Sec. IIIII this transition is first order. Second, con- 
sider disruption of Neel order by exchange randomness 
at T = 0. Strong disorder leads to a larger modulation 
of exchange interactions than that produced by a uni- 
form lattice distortion. In turn, this results in ground 
states with spin configurations that are random rather 
than Neel ordered, and an energy that is minimised at 
A = 0. The transition from the Neel ordered state takes 
place at S ~ We show in Sec. HVI that it is probably 
continuous. Finally, consider the effect of non-zero tem- 
perature at strong disorder. We have studied the model 
without magnetoelastic coupling (k — > oo) elsewhere, 
finding a phase transition from a spin-glass ordered low 
temperature phase to a paramagnetic high temperature 
phase, with a transition temperature T g ~ S. 9 We include 
this phase boundary in Fig. [2j although we have not at- 
tempted to investigate the multicritical point at which 
the three phases meet. 



III. MAGNETICALLY-DRIVEN STRUCTURAL 
PHASE TRANSITION 

In this section we discuss the thermally driven transi- 
tion in a system without disorder, between a high temper- 
ature, cubic, paramagnetic phase and a low-temperature, 
tetragonal phase, which we will show is Neel ordered. We 
restrict our attention to the regime of weak magnetoelas- 
tic coupling, k 3> J -1 , in which behaviour depends only 
on the variable kT and is independent of J. Our main 
tool is the use of Monte Carlo simulations to determine 
the free energy of the spin system. 

Let F(A, T) be the free energy per spin at tempera- 
ture T, calculated from the Hamiltonian Hq + Ha- The 
combination 



r(A,T) 



Hd 
~N~ 



F(A,T) - F{0,T) 



(6) 



is the net free energy gain per spin arising from a lattice 
distortion. The physical value of A is the one that min- 
imises r(A, T), and the system is in a tetragonal phase 
if there is a value of A > for which T(A, T) < 0. 
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FIG. 2: Schematic phase diagram in the plane of scaled tem- 
perature kT and disorder strength nS, based on numerical 
simulations of J5} . T c denotes the temperature of a transition 
between a Neel ordered, tetragonal low temperature phase, 
and a paramagnetic, high temperature phase in the pure sys- 
tem. A similar transition occurs at zero temperature, at a 
disorder strength S c . Solid lines represent phase boundaries 
that have been probed in simulations. Dashed lines are guides 
to the eye. 



The contribution of the spin degrees of freedom to this 
free energy for T, A <C J can be written in the scaling 
form 



F(A,T) - F(0,T) = A/ S (A/T) 



(7) 



since in this regime A sets the only energy scale for 
the spin system. The function f s (x) plays an important 
role in the following, and it is useful to discuss its low 
and high temperature asymptotics. At low temperature 
(i>1) the magnetic free energy difference is dominated 
by the energy gain 2A per spin in the Neel state, and 
so lim^^oo f s (x) — —2. Conversely, a high-temperature 
expansion (covering the regime 1 -C kT <C k J) gives for 
small x 



NAf s (A/T) = (H A )o - ^[{(Ha) 2 ),-{Ha)1\ 
+ C(A 3 /T 2 ), 



(8) 



where (. . .)o denotes an average with the Hamiltonian Ho 
in the limit T <C J. This average yields (Ha)o = and 

(Ha)q = bNA 2 with b > 0, and so f s (x) bx/2 for 

x -> 0. 

Now Eq. ([6]) can be written in terms of f s (x) as 



r(A,T) 



A 2 + A/ S (A/T) 

kT 2 

—x +xf s {x) 



(9) 



The form of the function / s (a;) determines the critical 
point kT c and the order of the transition. The equation 
kTx/2 + f s {x) = has no solution for positive x in the 
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cubic phase (kT > nT c ) but has a solution x c > in the 
tetragonal phase. If x c approaches zero as kT approaches 
kT c from below, the transition is continuous, while if x c 
remains finite in this limit, the transition is first order. 
Thus, the transition is continuous if f s (x) is convex for 
< x < oo; otherwise it is first order. 



A. Numerical Approach 

To find the magnetic free energy difference between the 
cubic and tetragonal phases, and specifically the function 
f s (x), we use Monte Carlo simulations of a spin system 
with the Hamiltonian Ho + H.A ■ The system simulated 
had the shape of a rhombohedron, with edges parallel to 
the primitive basis vectors of the lattice and with peri- 
odic boundary conditions. We focus on the temperature 
range T -C J for fixed A C J and ensure equilibration 
at low temperature by employing parallel tempering. 10 
Representative simulation parameters are A/J = ICP 3 
and T/J > 1CT 4 , with system sizes 32 < N < 2048 and 
measurement run lengths of up to 2 x 10 5 parallel temper- 
ing steps, preceded by equilibration for the same number 
of steps. Within the parallel tempering algorithm spin 
configurations are exchanged amongst a range of possible 
temperatures: we checked in particular that each config- 
uration explored all temperatures across an interval that 
spanned the first-order transition temperature. The data 
presented here are predominantly for systems of 500 and 
864 spins, but larger system sizes were also simulated to 
ensure that finite size effects are not significant. 

Determination of the free energy requires calculation 
of the entropy S(T). We do this by integration of the 
heat capacity C V (T), using 



S(T ) - S(T) = [ " 
Jt 



CAT) 
T 



dT 



(10) 



and choosing the reference temperature To ^ A so that 
S (To) is independent of A and cancels from the free en- 
ergy difference F(A,T) — F(0,T). We note that since 
C V (T) is constant at low temperature in our classical 
model, S(T) is logarithmically divergent as T — > 0, but 
the contribution — TS(T) to the free energy has no di- 
vergence. For fixed A > the spin system has a Neel or- 
dering transition, which is first order in the regime of in- 
terest, A <C J- Numerical integration of C V (T)/T across 
this transition is difficult, because the latent heat of the 
transition appears in a finite system as a narrow spike 
in C V (T). We circumvent this difficulty by using two in- 
dependent approaches to find the entropy change at the 
transition, and checking between them for consistency. 
One approach is direct integration, using spline fits of 
C V (T) around the transition. The other is to use the 
change in internal energy at the transition to determine 
the latent heat and hence the entropy jump. 



B. Results 

Results from simulations for the scaled free energy dif- 
ference [-FX A, T) — F(0, T)]/A and internal energy E per 
spin as a function of temperature are shown in Fig. [3] 
We discuss first the behaviour of E. For a system with 
A = the ground state energy per spin is — J, and E 
increases linearly with T at low temperature, because 
C V (T) is constant. For A > the ground state energy 
per spin is reduced by 2A. For fixed A there is a Neel 
transition, which takes place at T/A w 3.1: at the tran- 
sition the internal energy per spin increases abruptly and 
approaches its value in a system with A = 0. Consider 
next the the free energy difference. It is also small in 
the high temperature phase. In contrast to the internal 
energy difference, the free energy difference is of course 
continuous through the transition. It grows with decreas- 
ing temperature below the transition. 




FIG. 3: Scaled free energy difference (F(0) - F(A))/A vs 
T/A for A/J = 1CT 3 . Inset: internal energy per spin on the 
same temperature scale, for A/J = 10~ 3 (o) and A = (•). 
All data for N = 256 spins. 



The free energy difference is rcplottcd as a function 
of the variable x = A/T in Fig. [H together with the 
straight line through the origin that is tangent to the 
function f s (x). From the gradient of this line we identify 
the transition temperature T c in the system with mag- 
netoelastic coupling, obtaining nT c = 3.34. We also see 
that the transition is strongly first order, since the so- 
lution x c ~ 0.60 to the equation f 8 (x) — —kT c x/2 is 
finite. Note that the transition for the full model with 
Hamiltonian TL, in which the value of A varies with T to 
minimise the free energy, preempts the one in a system 
with fixed A. The latter transition is responsible for the 
kink in f s (x), visible at x » 0.3 in Fig. 2] Of course, it 
is the extra elastic energy accompanying a lattice distor- 
tion that destabilises the Neel ordered state at a lower 
temperature in the full model than in the one with A 
fixed. 

To demonstrate that the tetragonal phase is indeed 
Neel ordered, we show in Fig.[5]the variation of the sub- 
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1 2 3 4 5 

x=A/l 

FIG. 4: f„(x) vs x = A/T for fixed A/J = 10~ 3 (•). The 
straight line is tangent to f a (x) and has gradient —1.67. 



lattice magnetisation with temperature in a system with 
a fixed value of A. At the critical point of the full model, 
T/A m 1.7 and the sublattice magnetisation is approxi- 
mately 90% of its saturation value. We discuss the tran- 
sition at fixed A in detail elsewhere*AA 

In summary, the simple approach we have presented, 
which combines a microscopic treatment of the statistical 
mechanics of the spin degrees of freedom with mean field 
theory for the lattice distortion, yields a strongly first- 
order transition between cubic and tetragonal phases, 
with highly developed Neel order at all temperatures 
within the tetragonal phase. A first-order transition is 
also predicted using Landau theory, but the mechanism 
in that case is different: it arises because symmetry per- 
mits a term in the elastic energy which is cubic in strain. 
In our approach such a term would appear as a contribu- 
tion to Ti. e \ cubic in A. Within Landau theory, first order 
transitions can arise in the absence of cubic invariants, if 
the expansion of the free energy in powers of the order 
parameter has a quartic term with a negative coefficient. 
This appears to be the mechanism that leads to the first 
order transition we find. 



IV. INFLUENCE OF QUENCHED DISORDER 

As we have seen, a structural phase transition is a 
generic feature of the model we have studied, irrespec- 
tive of the strength of magnetoelastic coupling. The 
driving mechanism is the linear gain in exchange energy 
with lattice distortion for some ground state spin config- 
urations of the high-symmetry structure. This magnetic 
Jahn- Teller mechanism is likely to apply to a variety of 
other classical models for geometrically frustrated mag- 
nets. Specifically, for it to be effective the set of ground 
states should include configurations in which S^ • Sj takes 
different values for different nearest neighbour spin pairs, 
so that the exchange energy of the state varies linearly 
with lattice distortion A. For example, this condition in- 
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FIG. 5: Sublattice magnetisation m vs T/A with A/J = 10" 3 
for N = 864 (•) and 4000 (o). Errors only shown for N = 4000 
for clarity. 

eludes Heisenberg model on the lattice appearing in the 
material SCGO^ but excludes the kagome Heisenberg 
antifcrromagnet, in which Si • Sj ; = — 1/2 for all nearest 
neighbours in all ground states. In many materials, how- 
ever, no structural transition is observed. Instead, ge- 
ometrically frustrated magnets often show spin freezing 
at low temperature. Such freezing suggests the impor- 
tance of residual disorder, even though samples in some 
cases appear to be almost disorder-free. In this section 
we examine how disorder may stabilise the cubic phase 
of our model. We present results obtained in two ways. 
First, we use Monte Carlo simulations to study a spin 
system with fixed A > at various disorder strengths 5 
as a function of temperature. We find that strong dis- 
order suppresses the Neel transition that is observed at 
5 = 0. Second, we examine behaviour at T = as a 
function of disorder strength, treating the dependence of 
ground state energy on 5 in the same way as we did the 
dependence of free energy on T in Sec. IIIII 



A. Behaviour for T > at fixed A 

The introduction of random exchange interactions 
complicates Monte Carlo simulations. Longer equili- 
bration times are required (up to 10 6 parallel temper- 
ing steps) and we must average over disorder configura- 
tions. This limits the system sizes that can be simulated. 
The results we present were obtained for systems of size 
N = 864, averaging over 100 realisations for each disor- 
der strength, with A/J = 10 -2 . 

We show in Fig. [5] the heat capacity per spin as a 
function of temperature for various disorder strengths, 
expressed in terms of A/ 8 since we expect this ratio to 
control behaviour in the regime A, 5, T <C J. We find 
that the sharp peak in heat capacity, associated at 5 = 
with a discontinuous Neel ordering transition, occurs at 
lower temperature with increasing disorder but persists 
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FIG. 6: Heat capacity per spin vs T/J for A/S — 
10(D), 10), and 0.5(a) in a system of 864 spins at A/J = 
10 -2 . Errors shown for system with largest error. Inset: heat 
capacity per spin vs T/J for 500 (•) and 864 (o) spins at 
A/S — 10 _1 and A/J = 10 -2 . Errors indicated by scatter of 
points. 




FIG. 7: f E (y) vs y for three values of A/J: 10 _1 (o), 10" 2 
(o) and 5 x 10 -3 (□). The straight line is tangent to /e(j/) at 
y — and has gradient — kS c /2 — —1.6 The inset shows the 
same plot but for a larger range of y. Errors are less than the 
symbol sizes. 



for A/S > 1. At larger disorder strength, A/S < 1, it is 
rounded, being smooth and independent of system size 
at A/S = 1CP 1 (see inset to Fig. [6J. The dominant errors 
are from disorder averaging. We conclude that the Neel 
transition is suppressed by disorder. It is presumably re- 
placed at strong disorder by a spin glass transition, but 
because of the difficulties attached to Monte Carlo simu- 
lations in this regime, we do not investigate this further. 



B. Behaviour at T = 

To make additional progress in quantifying the effect 
of disorder on the structural and magnetic phase tran- 
sitions, we focus on behaviour at T=0 as a function of 
S. Let E(A, S) be the disorder-averaged ground state 
energy per spin for the spin system with Hamiltonian 
Wo + Ha + T~idis- In analogy with Eq. ©, we define the 
energy combination 



r E (A,J) = ^i 



E(A,5)-E(0,S). 



(11) 



To determine the ground state lattice symmetry for a 
given disorder strength S, the value of A should be chosen 
to minimise Te(A, <5). In the cubic phase the minimum 
is at A = 0, while in the tetragonal phase it is at A > 0. 

We expect the spin contribution to this energy differ- 
ence in the limit A, S <C J to have the scaling form 



E(A,S)-E(0,S)=Af E (A/S) 



(12) 



with the asymptotic behaviour for /e(j/) given at large 
y by lim^oo fs(y) = -2, and at small y by f^{y) oc y. 
Using this scaling form 



Te(A, 5) = A 



y2/ + /E(y) 



(13) 



Thus the state of the system is controlled by the dimen- 
sionless coupling constant nS. As for the thermally driven 
transition without disorder, there is a critical coupling 
kS c . If kS > kS c the equation nSy/2 + /e(?/) = has 
no solution for y > and the system is in the cubic 
phase. Alternatively, if kS < n5 c there is a solution at 
y c > and the system is in the tetragonal phase. As for 
the thermally driven transition, the disorder-driven one 
is continuous if /e(^) is convex for < x < oo and first 
order otherwise. 

We show in Fig. [7] the function fs(y) determined nu- 
merically for a system of size N = 864 spins, averag- 
ing over 50 disorder realisations. Low energy configura- 
tions are obtained by quenching low temperature states 
taken from Monte Carlo simulations using parallel tem- 
pering. The agreement between data for three values of 
A/J demonstrates that the calculations are within the 
scaling regime A, 5 <C J. In contrast to behaviour at 
the thermally driven transition, it appears from the form 
of /e(2/) that y c will approach zero smoothly as kS ap- 
proaches kS c from below. Hence the disorder driven, zero 
temperature transition is probably continuous. We find 
kS c = 3.2. 

As a supplement to this study of the structural tran- 
sition, we also examine the variation of sublattice mag- 
netisation m with disorder strength at fixed A in ground 
states. The sublattice magnetisation shows significant 
finite-size effects in our range of system sizes, and so we 
extrapolate to the limiting value TOoo from results for 
N = 256, 500 and 864. We show the variation of rriao 
with y in Fig. [51 and the extrapolation in the inset. Sub- 
lattice magnetisation is not much affected by disorder 
for y > 1, and decreases to zero for y w 0.17. If it is 
indeed the case that /e(2/) is convex, so that the struc- 
tural transition is continuous, and if it is also the case 
that the magnetisation is zero below a finite value of y, 
then the disorder driven, zero temperature transition is a 
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FIG. 8: Staggered magnetisation per spin in the infinite sys- 
tem vs y" 1 on a logarithmic scale. Inset: illustration of the 
extrapolation of sublattice magnetisation m vs N^ 1 to infinite 
system size. 



two-stage one, and there exists a tetragonal phase with- 
out Neel order at intermediate disorder strength. A more 
detailed investigation of this point would be demanding, 
and we do not pursue it. 



C. Paramagnetic-Spin Glass Phase Transition 

We have not investigated the regime k5 3> 1 as part 
of the present work, but we draw attention to earlier 
studies of the effect of weak disorder (8 <C J) at A = 
(in effect, k — > oo)&12- A continuous phase transition 
between the paramagnetic and spin glass phases occurs 
at a temperature proportional to 5. 9 



V. DISCUSSION 

The model we have discussed provides a caricature for 
the chromium spinel oxides ACr204, where A is Mg, 14 
Zn£ Cd,ii or Hg^ These 5 = 3/2 pyrochlore antifer- 
romagnets undergo structural phase transitions accom- 
panied by Neel ordering, with transition temperatures 
that are small compared to the magnitude of their Curie- 
Weiss constants. The structural distortions and mag- 
netic ordering patterns are different in each compound 
and are considerably more complex than the one we have 
considered; 1 -' 1 - While we are therefore not able to make 
a detailed comparison with experiment, some qualitative 



points deserve emphasis. These concern the mechanism 
driving the transition, the order of the transition, and 
the influence of non-magnetic impurities. 

Our calculations demonstrate that, at least in the clas- 
sical limit and for the uniform distortion mode we con- 
sider, the mechanism driving the structural transition is 
the energy gain from Neel order, rather than spin Peierls 
order. It is possible in principle that other distortion 
modes would allow a phase at intermediate temperature 
that has spin Peierls but no Neel order, but the very 
large discontinuity we observe in the Neel order parame- 
ter suggests to us that this is unlikely. It is also possible 
that quantum fluctuations favour spin Peierls order, and 
it would be interesting to attempt in future research to 
quantify this. 

In addition, our results show that the phase transition 
can be first order for a reason unconnected to the symme- 
try analysis of Landau theory. This is significant because 
within Landau theory some structural distortion patterns 
are expected to lead to a first order transition while oth- 
ers should not. In particular, within Landau theory a 
transition involving a tetragonal distortion that is pre- 
dominantly an odd mode (affecting the two tetrahedra 
in a unit cell oppositely) should be continuous, while a 
transition involving an even mode should be first order 2& 
Our simulations demonstrate that a first order transition 
may result from a microscopic treatment, independently 
of expectations based on symmetry. 

Finally, it is interesting to compare our study of the 
influence of exchange randomness on the transition with 
observations of the effect of non-magnetic impurities in 
the material Zni-zCd^C^O^ 17 In this system a sur- 
prisingly small level of substitution (x = 0.03) at non- 
magnetic sites is sufficient to suppress Neel order and the 
structural phase transition. The origin of this behaviour 
is believed to be the larger ionic radius of Cd 2+ compared 
to Zn 2+ : as a consequence, substitution leads to random 
strains, which in the presence of magnetoelastic coupling 
generate exchange randomness. The results we present 
in Sec. IIVI show this mechanism in action, and are qual- 
itatively consistent with the phase diagram observed in 
Zni_ a; Cd x Cr204 as a function of temperature and dis- 
order strength, with x playing the role of our variable 
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